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1. Introduction 

The concept of transverse-momentum-dependent (TMD), or fcj^-dependent, or "un- 
integrated" , parton distributions is frequently encountered in the QCD literature. 
This is so particularly in the small- a; domain where they play a prominent role. 
These distributions are to be contrasted with the integrated parton distributions 
(parton distribution functions, pdfs) which play a fundamental role in global QCD 
fits. 

For the fc_L-dependent distributions there is no unique definition, and as a conse- 
quence it is for the phenomenology of these quantities very important to relate the 
different definitions. We shall therefore here analyze some of the most commonly 
found definitions in the literature. 

Before going any further, we should mention that there exist two different and 
distinct sets of users of these distributions, each with corresponding domains of 
application: 

• One area we might characterize as "hard scattering factorization" , where the pri- 
mary object of attention is a hard scattering factorization property. An example 
is the Drell-Yan process at low transverse momentum. 

• The other area is that of small x, DIS with x -C 1 being a typical example. 
Here the emphasis starts with the BFKL formalism (or a generalization of it) for 
scattering processes in the Regge region. 

We are here primarily interested in the small- a; domain. The theory of the "hard 
scattering" domain is in great detail explained in p] (see also [2] for a short re- 
view) , and a compact presentation is also given in the talk by John Collins [3] . We 
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also note that the terminology "TMD" distribution is usually used in the "hard 
scattering factorization" while traditionally in the small- a; domain one speaks of 
"/c^-dependent" , or more commonly of "unintegrated" parton distributions. 

2. The parton model as background 

In the parton model, the concept of a parton distribution quantifies the intuitive 
expectation of a number density of partons of given flavor in the target hadron. A 
mathematically exact definition can be given in the light-front quantization as 

f . /h(xkl ) y 1 <^%U°fc.«,ji^) () 

/jMs>*-U-Z, 2a . (27r)a (p,h\P,h) ■ () 

Here j and h label parton and hadron flavor, a is a parton helicity index, P, h) is 
the target state of momentunE) P, and a and are parton light-front annihilation 
and creation operators. We use light-front coordinates, the hadron target has zero 
transverse momentum, and x = k + / P + . Integrating Eq. ([IJ over all k± one obtains 
the integrated distribution 

fj/h{x)= d 2 kj_f j/h (x,k±). (2) 

Jail k± 

The exactness of the relation, however, depends on the specific assumptions of the 
parton model which no longer hold in full QCD. 

A superficial glance over the relevant literature reveals that the number density 
interpretation of the TMD, fc^-dependent, and unintegrated parton distributions is 
taken rather literally also in QCD. The actual technical definitions in the different 
cases, however, may or may not conform to the idea of a number density. We 
also note that the parton model relation between the integrated and the TMD 
distributions in ([2]) is taken to be true in QCD as well (with the difference that the 
k± integral is performed only up to the hard scale Q). The question then is what 
exactly the justifications are for these assumptions. 

3. The gluon distribution at small- a; 
3.1. BFKL and the dipole picture 

The prototype of the QCD applications in the small- a; domain is the BFKL formal- 
ism [3H5]. Here the 7*7* scattering amplitude (for A + B — > A' + B') is in the Regge 
limit (s/t 3> 1 with t = —qj_) written as 

1mA(s,t) f d 2 k ± d 2 k' A 



k\ (k' x -q ± r 



I A (k±,q±)I B (k / ± ,q±)F(s, k±, k' ± ,q±), (3) 



a Strictly speaking the momentum states \P) do not belong to the ordinary Hilbert space of states 
since they are not normalizable. One can remedy this by replacing them with wave packets, and 
then taking appropriate limits as these packets approach the momentum states, but we need not 
worry about such level of rigor here. 
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where A and B denote the two virtual photons. The dimensionless objects I A and 
1b are called the impact factors, while F is the object satisfying the BFKL equation 
(Eq. (fTS)) below), and is commonly referred to as the "BFKL Greens' function". If 
by 4>a and 4>b one denotes the couplings of gluons to the external photons A and 
B (via quark boxes), then to lowest order in the singlet channel in Feynman gauge 
(two gluon exchange) the amplitude can be written in shorthand notation as 

~e lk I a Ib (0|T ( Pa ■ A) (p A ■ A) (p B ■ A) (p B -A)\0). (4) 

Here ~ e ik indicates the eikonal approximation in which the numerator of the Feyn- 
man gauge propagator g^ is replaced by p^p A /pA ■ Pb (PaPbIpa • Pb) for gluons 
coupling to A (B). Therefore to lowest order F is given by the vacuum expectation 
value of four off-shell gluons. To all orders, the gluon fields have to be summed into 
eikonal gauge links, i.e Wilson lines, which make F gauge invariant. 

In the space-time picture by Balitsky [7] , the qq pair emerging from the electro- 
magnetic currents J from one of the particles (say A) travels along a straight line 
trajectory and scatters off the gluon field created by B which has a delta function 
shape, i.e a Shockwave, in the high energy approximation. The propagation of the 
qq in the external gluon field created by B is represented by the Wilson lines 

U A (xx) = P exp f-ig J d\n A ■ A a (x± + Xn A )t F ^j , (5) 

where the color matrix t F is in the fundamental representation, and ha is a unit 
vector in the direction of p A . This picture also leads to the so-called dipole for- 
malism, as the qq pair emerging from the current J is in a color singlet state and 
thus can be seen as a "color dipole". The lowest order result in Eq. (Q} is in the 
formulation by Balitsky generalized to 

W g (M) = J d 2 k ± I A (k ± , qi _)(Tr{U A UA}(k^q^) (6) 

where 

(Tr{U A U A }(k ± ,q ± )) = J d 2 x ± e- lk ^^ (Tr{U A (x ± )U A (0±)}) . (7) 

The exact definition of the bracket (O) is given in [7] . Here it suffices to say that it 
involves a convolution of the Wilson lines with the currents Jb and Jb' ■ Notice that 
the Wilson lines here correspond to the all order summation of the factors pa • A in 
<U]), while the fields ps • A are kept implicitly in the definition of the bracket. When 
made explicit they should give rise to an additional pair of Wilson lines Ub in the 
direction of ps ■ 

In formula ^ the main object is the BFKL function F. From © we see that 
one has 

(Tr{U A U A }(k ± ,q x )) = J ] ^j^—^I B (k' x ,q ± )F(s,k x ,k , ± ,q x ). (8) 
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The expectation value of the Wilson line pair is essentially taken to be the "gluon 
distribution" (also called the "dipole gluon distribution") in the small- a; domain. 
While the formulas here are derived for 7*7* scattering, a similar identification 
of the "gluon distribution" is made for other processes as well, such as DIS off 
a proton or a nucleus, and also in proton-proton (pp), proton- nucleus (pA), and 
nucleus-nucleus (AA) collisions. For these latter processes we have not been able to 
find any proofs showing that this is indeed fully legitimate. The definition of the 
"dipole gluon distribution" is then generally taken to be 



where we leave the prefactor C unspecified since there does not seem to be a univer- 
sally accepted value in the literatur^. We should also mention that one needs to 
insert transverse gauge links at infinity to make the operators in ([6]) and ((9]) exactly 
gauge invariant. 

We notice that the structure of (j9]) does not really conform to the parton model 
definition (p} of a number density. Therefore, despite its name, this object does not 
correspond to the direct generalization of the parton model definition (fTJ. Moreover, 
the Wilson lines in the definition are always take in the fundamental representation 
while the generalization of the parton model result to full QCD naturally leads 
to Wilson lines in the adjoint representation for the gluon distribution [T]. It is, 
however, still common to assert the relation © to the integrated pdf, by writing 



By G we here denote the gluon pdf. For an explanation of the meaning of the 
parameter x in J-dip, see below. 

A reason for asserting this result is that the leading twist approximation of 
the BFKL evolution, which Jdip satisfies, reproduces for G, via (fTQ|) . the DGLAP 
evolution in the approximation where rif = and where only the 1/z term in the 
splitting function P gg is kept. 

Despite this agreement in the simplified limit, it is, however, a bad practice 
to hold on to the relation (fT0| for several reasons. One reason is that the partial 
agreement of the evolution equations does not generally imply that the relation 
(|10[) is a priori consistent with the respective operator definitions. The integrated 
pdf, G, which follows from the collinear factorization approach has a clear operator 
definition (see for example PQ) that cannot simply be reduced to the integral of 
(j9|) (or some variant of it). Moreover, it is generally for phenomenology important 
to include the non-singular parts of the DGLAP splitting function, even at rather 
small-x. Therefore the simplified limit mentioned above is not in practice very useful. 

b In fact we have find several incompatible values of C used in different studies of the same processes. 
Moreover, even the form of varies from case to case so it is hard to pin down any exact universal 
definition. 






(10) 
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There are also conceptual problems. It is known that the use of light-like Wilson 
lines in the operator definitions for the TMD distributions gives rise to rapidity 
divergences which have to be somehow cut-off. A convenient way to cut the diver- 
gences is to take the Wilson lines to be non-light-like [U[7], i.e. to let them have 
finite rapidity. It is actually this cut-off which gives rise to the evolution in rapidity, 
just as the cut-off /i in momenta gives rise to the standard DGLAP evolution with 
respect to u. If we denote the cut-off in rapidity by £, then J-"di p depends on k± 
and £, J~dip = J 7 di P {k±; £). It is important to realize that £ is conceptually different 
than the variable x appearing in fj/h( x , k±) in the parton model definition (fTJ). In 
that case x is an "intrinsic" or kinematical variable, and is literally the light-cone 
momentum fraction x = k + / ' P + of the parton which participates in the hard scat- 
tering. On the other hand £ is an arbitrary cut-off variable which determines the 
total rapidity range for the soft gluons associated with J-dip, and has no counterpart 
in the parton model. Now, in the small-x literature one always sets Q — x (we remind 
that in DIS x — Q 2 /2q ■ P) . This is the reason behind the notation J-"di P (a;, k±). 

The variable x in the integrated distribution is not really related to the rapidity 
cut-off in J-"di p . It is in the integrated distribution a kinematical variable which has 
the meaning of the total light-cone momenta k + = x P + exchanged in the i-channel 
between the hard scattering and the target hadron. On the other hand in J-"di p it 
arises as the rapidity cut-off which can be implemented by taking the Wilson lines 
U to be non light-like. In fact, the corresponding kinematical variable x has in the 
derivation of (J3j) , ((6]) and (|9]) already been set to (that is, the direct analogue of 
x in ([1]) has already been set to in ©). A more correct notation for J^ip would 
therefore be J-di V {x = 0,£;_l;£ = x). The integrated distribution G does not have 
any £ dependence because it does not contain any rapidity divergences. We therefore 
see that one generally has to be more careful when relating G and J"di p . 



3.2. The "Weizsacker-Williams" distribution 

There is also a different type of gluon distribution found in the small- a: literature 
which is meant to literally be the analogue of the parton model definition (TTJ) . This 
distribution has been dubbed the "Weizsacker-Williams" (WW) gluon distribution, 
and is in a sense somewhat more closely related to the gluon distribution in the 
hard scattering factorization approach. A definition is provided in [5] which starts 
from the relation 

^ = (a i a ^x+,k)a i a {x + ,k)) = -^(Al(x+,k)Al(x + ,-k)) (11) 

where a and are the light-front gluon field annihilation and creation operators. 
Noticing that in light-cone gauge, A + = 0, one has F+ l = ik + A l a , it is seen that 
this definition thus corresponds to the expectation value (F+ Z F+ Z ). Indeed this is 
the direct generalization of the parton model definition ([1]) . The result in a generic 
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gauge is then taken to be [5] 

JW(z,fcx;C) = ^3 J dx~dy~J d 2 x x d 2 y x e ixP+ ^-^- ik ^- y ^ 

(F+\x)U ab {x,y)F+\y)) (12) 

Here U(x, y) denotes a Wilson line from y = (0 + ,y~,y±^) to x = (0 + , x~ , x±) in the 
adjoint representation. 

It should be noted that the definition (fT2j) is classical in the sense that all gluon 
fields in the expectation value are classical color fields, and the averaging (0}w p+ 
is performed using the classical distribution function W of CGC. All the effects of 
quantum fields has been inserted into the evolution of W with respect to the cut-off 
parameter £ (indicated by the notation W^p+). The role of £ is here identical to 
the cut-off discussed above. In fact, in the CGC, (P + acts a sharp cut-off so that 
gluons with k + < (P + are removed completely. Note that we again keep distinct 
the variables x and £. 

Notice that the mere requirement of gauge invariance does not fix the paths of 
the gauge links in (fT2")l exactly. Thus if the general form of the result is only dictated 
by gauge invariance, then the paths of the Wilson lines remain ambiguous. In the 
hard scattering factorization, the Wilson lines of the TMD distribution are fixed 
by whatever is needed in order to obtain factorization. In the CGC, however, the 
distribution (|12[) is not related to any factorization theorem. Its existence is simply 
conjectured from the intuitive expectation of the gluon distribution corresponding 
to a number density. It is furthermore not clear how (|12p arises from Feynman graph 
calculations in the small- a; limit since there does not exist any direct derivation of 
it. 

As the WW distribution is in structure more similar to the TMD distributions 
in the hard scattering factorization approach, it would be interesting to actually 
find a derivation of it from a factorization theorem in the small- a; limit. For some 
related work see also [9]. 



3.3. The CCH approach 

The standard reference to /c^-factorization in small- a; physics is the work by Catani, 
Ciafaloni and Hautmann (CCH) [TO]. One studies here the production of heavy qq 
pairs in photoproduction, DIS, and in hadron-hadron collisions, and the main goal 
is to formulate a TMD factorization formula at small- a; which at the same time can 
in the collinear limit be related to the standard collinear factorization formula. 

The factorization formula, based on the properties of the light-cone gauge, is 
written as (in case of photoproduction) 

(7 ~ 19 = 2~KP J ~z~ J d2k ^Pl z > k ^ M ^ T ^ fc ^' ( 13 ) 
where p = AM 2 / s, and M is the invariant mass of the heavy quark. This formula 
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is somewhat different than the BFKL formula ([3]), and is more close in spirit to 
the hard factorization approach since a here plays the role of a hard scattering 
coefficient. The object referred to as the "unintegrated gluon distribution" can in 
the light-cone gauge used in [TU] be written as [TT] 

n**±) = J dX { ]^t X 4p+ e^*--*-°-(P\l¥(0 + ,x-,x ± )F+ i (p)\P)- (14) 

This definition exactly corresponds to the parton model one in (jTJ) (to xf g /h), and 
is also the analogue of the WW distribution in the CGC. If taken literally, however, 
it leads to rapidity divergences (in this case due to the singular denominator of the 
light-cone gauge propagator), and consequently it must be modified. The point is, 
however, that the definition is never actually used in |10) . It is instead stated that 
J 7 (z,k±) in (|13p is "defined" by the BFKL equation 

T(x,k ± ,Q 2 ) = -6(kl-Q 2 ) + a s f ^[ ^(^( x / z ,k ± + q ± ,Q 2 )- 

9(kl-ql)F(x/z,k ± ,Q 2 )^j, (15) 

where a s = a s N c /n. This switch implies that the rapidity divergence is cut off since 
there is an implicit cut in the BFKL formalism. Effectively one introduces a cut-off 
£, and then sets C = z in (JT3J) . 

This switch implies a conceptual change which is rather important to under- 
stand. As we have noticed, the would be definition (fT4|) in light-cone gauge cor- 
responds to the WW definition. However, once it is instead declared that T is 
"defined" via the BFKL equation (1X5]) . one comes back to ([3]) and essentially to the 
dipole definition ©, and the latter does not conform to the parton model idea of a 
parton distribution. There is therefore an implicit change in the meaning of T . Ac- 
tually, T is now instead replaced by F in (J3J, that is the "BFKL Green's function" 
which itself is not a "gluon distribution". It is in fact stated in [TU] that for any 
realistic calculation one should rather provide a non-perturbative Qq -distribution 
of the gluon in the hadron to be convoluted with ^(x, fcx; Qq)- What this simply 
means is that one needs to perform a convolution just as in the right hand side of 
Eq. ([8} (Qo here corresponds to k' ± in JS])) . 

It is clear, however, that this cannot be directly compatible with the definition 
([141) that follows from (TT31) , because the operator definition of ([9]) cannot be reduced 
to ([14]) in light-cone gauge. There is moreover another problem. In the set up of [10] 
which gives (|13|) . the photon actually scatters off a single gluon with momentum P. 
The state \P) in (|9]) therefore strictly speaking corresponds to an on-shell gluon with 
momentum P, and not to a hadrorH The idea is based on the so-called "factorization 
of mass singularities" which has been used when dealing with collinear factorization 

c The cross section in Eq. (1 1 3 I t is therefore really a partonic cross section which is indicated by 
the subscript jg. 
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(see for example [H]). In this approach it is first asserted that a hadronic structure 
function, \Y hadror \ is a convolution of the corresponding partonic structure function 
W parton and & so _ ca u e( } "bare parton density" G bare , 

W hadron (q, P) = W parton {q, P = £P) ® 4 G bare (0- (16) 

The convolution here in the longitudinal variable £ is the same as the z convolu- 
tion in (TT51) . Factorization is then understood as the procedure of extracting out a 
divergent factor D from \yP arton to "define" a "renormalized parton distribution" , 

w hadron = (fr g> £)) g, Qbare = - g, (£, g, G bar e ) = - ^ gre^ ^ 

For the conceptual problems of this approach in collinear factorization we refer the 
reader to pQ. In this case the procedure is thus extended to TMD factorization. 
However, the parameter Qq in (|15l) implies that the incoming gluon P is no longer 
on shell, and the type of convolution ([8]) which is supposed to give a definition of 
the "unintegrated gluon distribution" clearly is rather different than the procedure 
in Eqs. (fl6|) and (fl7|) . It is also clear that F in (|3]) is different than T defined in 

The question here should be: Can we in the small- a; region indeed formulate a 
formula like in (fl3|) and then show that it leads to a precise operator definition of the 
TMD gluon distribution which is related to the WW distribution, with all the sub- 
tleties regarding the rapidity divergences taken properly into account? Additionally, 
can we then properly relate this distribution to the integrated gluon distribution? 
Moreover, what evolution equation in rapidity will this distribution satisfy? In the 
TMD approach the rapidity evolution is given by the so-called CSS evolution pQ, 
while in the formalisms studied here by the BFKL equation. One should investigate 
the exact connection between the two. 

4. Speakable and unspeakable 

As we have seen above, the idea and concept of the "unintegrated gluon distribution" 
in small- a: QCD is still intuitively based very much on the parton model definition 
of a number density, but we have also seen that the "dipole gluon distribution" does 
not conform to this idea, even if it is intuitively thought of being so (and repeatedly 
said being so). The parton model obviously provides a very useful intuitive guidance 
in full QCD, but unnecessary confusion, and even wrong results, can easily arise if 
one is not careful. 

In the historical development of QCD it was of course the approximative Bjorken 
scaling observed in the SLAC-MIT data, together with the successes of the quark- 
gluon model, which eventually led to the identification of partons with quarks and 
gluons, thus providing the latter an ontological commitment as real entities (rather 
than simply being abstract mathematical objects). 

The scaling violations of QCD modify the naive parton model and the definition 
of the parton distribution functions which in the parton model are strict number 
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densities according to ((!]). The parton model is, however, not completely thrown 
away, but the full QCD result rather corresponds to an improvement, without too 
dramatic or violent differences. Very important in this is, however, the existence of 
factorization theorems in QCD (for a number of elementary processes). For assume 
that factorization did not hold (its existence requires after all non-trivial proofs) . In 
that case the concept of parton distributions would not be useful at all, and it would 
not be possible to gain any real knowledge of the underlying entities, the partons, 
using them. In particular since partons are not directly observable due to confine- 
ment, we would be in a rather difficult situation in the detailed exploration of the 
inner structure of hadrons. In that case it would obviously be very questionable to 
commit any ontology to the partons, for example speaking of a "number density" , or 
a "fcj^-distribution" or the "probability of finding partons with momentum fraction 
between xP + and (x + dx)P + " and so on. Even when factorization does hold, how- 
ever, the strict number density interpretation is lost due to the UV renormalization 
in QCD, and the legitimacy of the statements made upon the parton distributions 
must be solely based on the mathematical rules of QCD. 

Quantum mechanics teaches us to speak of only what we can measure. In accel- 
erator experiments, what we measure are a bunch of hadrons and leptons in tracking 
chambers and calorimeters. To gain any knowledge of the partons from the observed 
hadrons we employ factorization and parton distributions. The former is, however, 
the crucial ingredient in this procedure, which is the moral of the story. In short we 
may say that parton distributions have no ontological priority over the factorization 
theorems. 

5. Summary and Outlook 

The problem of factorization is more intriguing and complex in the case of TMD 
factorization which contains more information on the underlying entities, but which 
also poses more difficulties. The fact that factorization appears to be broken in pp 
collisions [T3l[T4] (and mostly likely so in nucleus collisions as well) calls for caution 
in the applications of naive factorization. As the LHC is a pp (and AA) machine at 
very high energies, it is particularly crucial to understand the issue in the small-x 
domain. A preliminary overview of the literature shows a lack of proofs, and also 
some confusion as to the meaning of the "unintegrated gluon distribution" . 

There has recently been many applications of the "dipole gluon distribution" 
in single inclusive hadron production (or rather gluon production) , see for example 
[rSHTT] . The recent multiplicity measurements of ALICE [18] has shown that the 
central multiplicity grows faster with s in A A collisions than in pp collisions. The 
question is of course why this is the case. We must then be very cautious, however, 
since applications of small-x physics all assume the existence of factorization in these 
processes, and this is highly questionable. Moreover it is not at all clear that it is 
the "dipole gluon distribution" which is the relevant object to use in these cases. 
Given the enormous complexity especially of the AA collisions, it is no wonder that 
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the naive application of "dipole gluon distribution" does not work. 

We have undertaken a comprehensive study in the hope of provoking further 
work to address the issues mentioned here, so that the important question regarding 
factorization and the TMD parton distributions may be hopefully clarified. 
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